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ABSTRACT 


A  1 /24th  scale  model  of  the  150  FT.  WYTM  was  equipped  with  an 
air  bubbler  system.  The  bubbler  system  was  operated  at  various  flow 
rates  from  no  flow  to  twice  the  nominal  system  capacity.  Tests  were 
performed  in  both  solid  ice  and  brash  to  determine  the  effectiveness 
of  the  bubbler  system  in  reducing  the  power  required  to  transit  ice 
fields.  In  all  cases,  a  flow  was  found  which  reduced  the  total  power 
from  that  required  for  a  ship  without  the  air  system. 


NOMENCLATURE 


Beam 

Coefficients  of  Regression  Equations 

Elastic  Modulus  of  the  Ice  Sheet 

Mass  Moment  of  Inertia 

Length 

Mass 

Power 

Volume  Flow 

Resistance  or  gas  constant 

Temperature 

Work 

Coefficient  of  Kinetic  Friction  (friction  factor) 
Acceleration  due  to  Gravity 
Ice  Thickness 

Ratio  of  specific  heats  for  a  Gas 
Characteristic  Length  of  the  Ice  Sheet 
Pressure 

Time,  thrust  deduction  fraction 
Speed 

Thrust  identity  wake  fraction 
Ship-Model  Geometric  Scale  Ratio 
Mass  Density  of  Ice 
Mass  Density  of  Water 
Flexural  Strength  of  Ice  Sheet 
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1 .  SUMMARY 


The  purpose  of  this  series  of  tests  was  to  determine  what 
airflow  to  the  air  bubbler  system  on  the  Coast  Guard  140-FT  WYTM  is 
most  effective  in  increasing  icebreaking  performance.  The  bubbler 
system  was  modeled  for  1.5  ft.  of  solid  ice  and  3.0  ft.  of  brash 
ice.  Three  speeds,  1,  3,  and  5  knots  full-scale,  plus  static  resis¬ 
tance  were  examined.  A  range  of  six  flow  rates  was  covered  from  the 
least  flow  required  to  decrease  resistance  to  the  point  where  in¬ 
creasing  flow  had  little  additional  effect.  The  most  effective  air¬ 
flow  is  that  representing  the  lowest  total  horsepower,  adding  air 
compressor  power  to  shaft  horsepower  at  each  of  the  three  speeds  and 
in  static  conditions.  The  following  Table  summarizes  test  results. 


BUBBLER  FLOW  FOR  LEAST  REQUIRED  HORSEPOWER 


Bubbler  Flow  Total 


Ice  Type 

Thickness 

Speed,  Kts 

SCFM 

Horsepower 

Solid  Ice 

1.7  ft. 

0 

5,000 

1 ,770** 

1.5  ft. 

1 

7,500 

2,050 

3 

6,500 

2,900 

5 

6,000 

3,450 

Brash  Ice 

3.5  ft. 

0 

4,000 

600 

6.0  ft.*** 

1 

13,000 

1,870 

3 

11,500 

2,120 

5 

8,500 

2,250 

The  performance  of  the  140-FT  WYTM  in  solid  ice  without  an 
air  bubbler  system  has  been  determined  in  earlier  tests  [1]*.  In 
these  earlier  tests,  the  change  in  resistance  with  varying  hull-ice 
friction  was  also  examined.  In  the  present  tests,  a  constant  hull- 
ice  friction  of  0.23  was  used. 

The  140-FT  WYTM  bubbler  system  is  designed  for  an  airflow  of 
7500  SCFM  [2].  Manifold  cross  sections  and  orifice  sizes  are  based 
on  this  flow.  Variations  in  manifold  pressure  with  flow  rate  are 
found  in  the  design  report  [2],  The  evaluation  above  is  based  on  a 
constant  piping  system. 

In  addition  to  the  collection  of  resistance  data,  this  pro¬ 
gram  also  included  high-speed  underwater  movies.  These  films  may  be 


♦References  listed  in  Section  6. 

♦*  Weak  Ice 

***  Extrapolated  from  test  conditions 
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used  to  evaluate  the  placement  of  nozzles  in  planning  for  possible 
different  arrangements  on  new  designs,  such  as  the  proposed  Great 
Lakes-Arctic  East  Icebreaker.  They  clearly  show  air  flow  patterns 
and  may  be  used  also  to  determine  if  air  will  be  injested  into  seachests. 
The  movies  indicate  that  additional  nozzles  further  forward  might  be 
helpful.  Seachests  located  near  the  keel  appear  unlikely  to  injest 


2.  MODEL  AND  TEST  DESCRIPTION 


2.1  Model 

A  wooden  model,  constructed  to  the  molded  lines  of  the  140 
FT  WYTM.was  used  in  these  tests.  This  same  model  was  used  for  the 
tests  conducted  earlier  [1].  Table  2.1  contains  the  particulars  of 
the  full-scale  tug  and  the  1 /24th  scale  model  while  Figure  2.1  is  an 
abbreviated  version  of  the  lines.  In  addition  to  the  rudder,  an 
appendage  representing  the  bubbler  duct  below  the  keel  forward  was 
added  to  the  model.  This  was  added  to  ensure  that  bubbler  flow  did 
not  cross  the  keel  since  orifice  placement  is  slightly  higher  than 
the  full-scale  nozzles.  Figure  2.2  shows  the  difference  between 
orifice  placement  on  the  model  and  on  the  full-scale  ship  for  the 
bow  manifolds. 


2.2  Bubbler  System 

Each  of  the  four  manifolds  on  the  ship  was  separately  modeled 
and  instrumented.  All  orifice  sizes  were  1/16  in.  (1.6  mm)  which 
represents  1.5  in.  full-scale.  Figure  2.3  is  a  schematic  of  the 
model  bubbler  system.  The  flowmeters  utilized  were  Brooks  Instru¬ 
ment  Model  1355  Sho-Rate  150's  on  each  manifold.  These  were  factory 
calibrated  tapered  tube  flowmeters.  Pressure  was  measured  using 
Dwyer  Instrument  Series  2000  Magnehelic  pressure  gages  immediately 
downstream  of  the  flowmeters.  Temperature  was  measured  in  the  mani¬ 
folds  utilizing  copper-constantan  thermocouples.  Manifolds  were 
numbered  as  shown  in  Figure  2.3  and  data  recorded  for  each  manifold 
separately. 


2.3  Arrangement 

Figure  2.4  shows  the  basin  layout  used  to  accomplish  the 
tests.  The  model  was  first  tested  on  the  East  Side  of  the  model 
basin  and  then  moved  to  the  West  Side.  The  airflow  and  velocity  con¬ 
ditions  tested  in  solid  ice  were  repeated  in  brash. 

2.4  Model  Scaling  Laws 

For  model  testing  to  be  correct,  the  model  must  be  both 
geometrically  similar  and  dynamically  similar  to  the  full-scale  pro¬ 
totype.  The  first  is  achieved  by  scaling  all  dimensions  by  the 
geometric  scale  factor  X.  The  second  condition  is  achieved  by  main¬ 
taining  the  ratio  of  significant  forces  the  same  for  both  the  model 
and  the  prototype. 

In  testing  a  ship  model  in  ice,  the  significant  forces  are 
gravity  forces,  dynamic  forces,  ice  forces,  and  friction  forces. 
Gravity  forces  will  scale  by  X3  since  the  density  of  water  is  the 


TABLE  2.1  SHIP  AND  MODEL  DATA 
for 

140  FT.  WYTM 
Model  No.  202C 

APPENDAGES  '•  Rudder  and  Bubbler  Duct  LINES  DWG  NO.  :140-WYTM-0101  -1 


DIMENSIONS 

ITEM 

SHIP 

MODEL 

Length  (LOA),  m 

42.67 

mm 

Length  (LWL),  m 

39.62 

■EH 

Length  (LBP),  m 

39.62 

wsm 

Beam,  Max  (at  LWL),  (3^) ,  m 

10.41 

0.434 

Beam,  Max,  m 

11.43 

0.476 

Draft  (at  test  WL),  m 

3.658 

0.152 

Trim  (test)  (+  AFT) ,  m 

0 

0 

Displacement,  Metric  Tons 

673 

0.0487 

Wetted  Surface  Area,  Sq.  m 

409 

0.710 

Distance  from  FT  to  3  ,  m 

21.79 

0.908 

Distance  from  to  LCF  (+  AFT),  m 

-1.393 

-0.058 

Distance  from  B  to  HI  (+  AFT),  m 

-1.981 

-0.083 

Distance  from  LCG  to  31  (+  AFT),  m 

-0.710 

-0.030 

Longitudinal  Metacenter  Above  Keel,  m 

5.15 

0.215 

COEFFICIENTS  AND  ANGLES  (AT  TEST  WL) 


Block  Coefficient 

0.437 

Length/Beam 

3.81 

Midship  Section  Coefficient 

.748 

Beam/Draft 

2.85 

Prismatic  Coefficient 

.589 

1.44 

Waterplane  Coefficient 

.715 

n2 

3.03 

Bow  Stem  Angle,  Deg. 

30° 

X 

24 

Half  Entrance  Angle,  Deg. 

27° 
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ABBREVIATED  LINES  AND  BOOT  PLAN 
0?  UNITED  STATES  COAST  GUARD  140- FOOT  WTTM 


FIGURE  2.1. 


Figure  2. 3.  Bubbier  System  Schematic. 
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Figure  2.4.  Basin  Layout 


same  for  both  the  model  and  the  full-scale  ship.  It  then  follows 
that  the  dynamic  forces  must  also  scale  by  X3,  which  is  achieved  by 
testing  the  model  at  speeds  corresponding  to  the  full-scale  speeds 
divided  by  /X.  This  results  in  the  Froude  number  v/JgL  being  the 
same  for  both  the  model  and  the  prototype. 

From  the  principle  that  all  forces  acting  on  the  model  must 
scale  by  X3,  the  scaling  laws  listed  in  Table  2.2  can  readily  be 
derived.  These  laws  dictate  that  the  model  ice  thickness  h,  flexural 
strength  erf,  and  elastic  modulus  E  must  be  reduced  from  the  appro¬ 
priate  full-scale  values  by  the  scale  factor  X,  and  the  density  of 
Pi  must  be  equal  to  the  full-scale  value.  Bubbler  scaling  laws  are 
discussed  in  Appendix  E. 

In  using  this  procedure,  the  viscous  forces  do  not  scale  pro¬ 
perly.  These  forces,  however,  will  be  negligibly  small  compared  to 
the  other  forces  involved. 

2.5  Model  Testing  Facility 

The  model  test  series  was  conducted  in  the  ARCTEC  Ice  Model 
Basin,  Columbia,  Maryland.  This  facility  consists  of  a  refrigerated 
model  towing  basin  30.5  metres  long,  3.7  metres  wide,  and  1.5  metres 
deep,  which  is  filled  with  a  saline  solution.  On  the  surface,  ice 
is  frozen  to  a  thickness  equal  to  the  geometric  scale  of  the  desired 
full-scale  ice  thickness.  By  controlling  the  water  salinity,  freezing 
rates,  and  temperatures,  model  ice  is  produced  with  properties  correctly 
scaled  for  model  testing,  with  the  exception  of  the  elastic  modulus. 

The  models  are  towed  through  the  ice  sheet  at  constant  speed 
and  the  resistance  is  measured.  The  carriage  drive  system  is  capable 
of  towing  at  any  desired  model  speed. 

2.6  Resistance  Test  Procedure 


The  two  test  tracks  were  arranged  side-by-side  in  order  to 
obtain  the  greatest  number  of  data  points  from  each  ice  sheet.  This 
procedure  has  been  shown  to  be  valid  providing  the  distance  between 
the  models  is  greater  than  six  characteristic  lengths  and  the  distance 
from  the  models  to  the  basin  walls  is  greater  than  three  character¬ 
istic  lengths  [1].  For  this  criterion,  the  characteristic  length  l 
of  the  ice  is  defined  by:  ° 


l 

a 


(2.1) 
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TABLE  2.2  SCALING  LAWS  FOR  MODELING* 


Variable 


Symbol 


*  Subscripts  fs  =  full-scale 
ms  *  model -scale 


Scaling  Law 


Length 

L 

Lfs 

3 

XLmc 

ms 

Beam 

B 

Bfs 

= 

xsmc 

ms 

Ice  Thickness 

h 

hfs 

= 

X*ms 

Ice  Flexural  Strength 

°r 

affs 

= 

Ice  Elastic  Modulus 

E 

Efs 

= 

AS 

ms 

Speed 

V 

ufs 

= 

Ay 

ms 

Time 

t 

*fs 

= 

Afc. 

ms 

Resistance 

R 

*fs 

= 

ms 

Mass 

M 

*Ts 

= 

a3a? 

ms 

Mass  Moment  of  Inertia 

I 

Jfs 

= 

A5I 

ms 

Mass  Density  of  Water 

PW 

pufs 

3 

pums 

Mass  Density  of  Ice 

p; 

pifs 

= 

pims 

Coefficient  of  Kinetic 
Friction 

f 

ffs 

= 

f 

Jms 

Volume  Flow 

Q 

«fs 

= 

X2,5«n 

Power 

p 

Pfs 

= 

X3‘5pn 

where 


E  -  elastic  modulus  of  the  ice  sheet 
h  =  thickness  of  the  ice  sheet 
p  =  mass  density  of  water 
g  =  acceleration  due  to  gravity 

The  ratio  d/l„  (where  d  is  the  distance  from  the  side  of  the  model  to 
the  basin  wall)  must  exceed  3.  Thus,  maximum  ice  thickness  for  dual 
model  tests  is: 


where : 

E  -  4  x  106N/m2  (Approximate  actual  basin  value,  not  correctly 

scaled) 

Pw  =  1005  kg/m3 

g  -  9.807  m 

d  =  543  mm  (see  Figure  2.4) 

All  tests  met  this  criterion.  By  using  approximately  two  model  - 
lengths  per  data  point,  the  models  could  be  towed  at  six  different 
velocities  in  each  ice  sheet,  thereby  collecting  twelve  data  points 
per  ice  sheet. 

For  each  test  the  speed  v  of  the  towing  carriage  and  the  re¬ 
sistance  R  of  each  model  were  recorded  on  an  oscillograph.  In  addition, 
the  flexural  ice  strength  Of  was  measured  in  three  locations  before 
and  after  each  track.  The  elastic  modulus  E  was  measured  at  three 
positions  prior  to  each  track, and  the  ice  thickness  was  measured  every 
metre  on  both  sides  of  the  broken  channel  following  the  test  run. 

The  throttling  valves  ahead  of  the  flowmeter  were  adjusted  until 
test  values  of  flow  existed.  Table  2.3  contains  a  list  of  test  values. 
After  two  ice  sheets,  the  lowest  two  values  were  added  and  the  high¬ 
est  two  dropped  when  there  appeared  to  be  little  change  with  the 
highest  flows.  Gas  temperature  and  pressure  were  measured  to  correct 
this  flow  to  standard  conditions  before  conducting  the  data  analysis. 

The  equation  used  to  correct  flowmeter  reading  was: 
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TABLE  2.3 

TARGET  AIR  FLOWS 

UTILIZED 

Name 

Full-Scale 

ft3/min 

Model -Scale 
ft3/h 

Bow  Man. 
ft3/h 

Amidships 

Man. 

ft3/h 

Qo 

-  NO  FLOW  -  - 

938 

20.0 

4.5 

5.5 

«>  A 

1,875 

39.9 

9.0 

10.9 

Qi 

3,750 

79.8 

18.1 

21.8 

*Qz 

7,500 

159.5 

36.1 

43.6 

Q3 

11,250 

239.2 

54.2 

65.4 

Qk 

15,000 

319.0 

72.2 

87.2 

*  Design 
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where : 


Qs  -  flow  at  standard  conditions 

Qn  =  flowmeter  volume  (instrument  error  removed) 

xl 

PN  =  nitrogen  pressure 

p  =  standard  pressure  =  14.7  psia  =  407  in.  H20 
Tg  -  standard  temperature  =  289°K  =  520°R 
r  -  nitrogen  temperature 

fi.  lL/£  \ 

rn  -  nitrogen  gas  constant  =  55.1 

i?air  =  constant  for  air  =  53.3 


All  data  is  tabulated  in  Appendix  A. 

The  resistance  of  the  model  was  measured  using  a  strain- 
gaged  force  block.  The  model  was  attached  to  the  force  block  in  such 
a  way  to  allow  pitching,  rolling,  and  heaving  motions.  The  models 
were  restrained  in  yaw  and  sway.  A  daily  calibration  of  the  force 
block  was  performed  in  order  to  ensure  accurate  measurements. 

The  model  speed  and  position  were  measured  simultaneously  by 
recording  the  passing  of  six  spokes  in  a  wheel  of  the  carriage  drive 
system.  Each  pulse  on  the  oscillograph  indicates  a  carriage  movement 
of  a  fixed  distance.  By  recording  the  distance  traveled  on  a  time- 
based  recorder,  the  velocity  can  be  calculated.  The  carriage  position 
in  relation  to  the  ice  sheet  was  determined  by  noting  the  starting 
position  of  the  models  and  then  counting  the  pulses. 

High-speed  underwater  movies  were  taken  at  a  frame  speed  of 
24A  frames  per  second.  When  projected  at  normal  speed,  the  motion 
of  the  model  is  viewed  in  full-scale  time.  Supplemental  surface 
footage  was  taken  with  a  normal  speed  camera.  Test  films  were  edited, 
spliced  together,  and  titled  to  produce  a  film  summary  of  the  program. 


Brash  Preparation  and  Monitorinc 


After  testing  in  level  ice,  the  ice  sheet  was  broken  up  into 
small,  random  size  pieces  of  which  50£  were  less  than  25  mm  on  the 
largest  dimension.  The  broken  ice  was  compacted  into  an  area  of  even 
thickness  of  30  to  45  mm  corresponding  to  an  average  full-scale  thick- 
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ness  of  3  feet.  Thickness  was  measured  using  the  gage  shown  in  Figure 
2.5.  This  modeled  unconsolidated  brash  Ice.  In  order  to  further 
characterize  the  model  brash  Ice,  two  additional  measurements  were 
made. 


The  first  was  made  with  a  cone  penetrometer.  This  consisted 
of  a  weighted  cone  made  of  transparent  plastic  with  a  40°  angle. 
Graduations  for  depth  of  penetration  were  marked  on  the  inside  of  the 
cone.  This  was  used  primarily  to  ensure  that  the  brash  ice  was  of 
uniform  compactness.  Figure  2.6  shows  the  cone  penetrometer. 

The  second  measurement  determined  the  mass  of  ice  per  unit 
volume.  This  was  made  using  a  wire  basket  (300  x  300  x  100  mm)  to 
lift  the  ice  out  of  the  water.  After  it  was  allowed  to  drain,  the 
weight  of  the  ice  in  the  basket  was  measured  using  a  shipping  scale. 


3.  RESULTS 


3.1  Resistance  Data 


All  of  the  resistance  data  collected  during  the  tests  is 
contained  in  Appendix  A.  There  are  three  sections  to  Appendix  A. 

The  first  section  contains  the  resistance  data  using  the  units  in 
which  it  was  measured  at  the  model  scale.  The  second  section  of 
Appendix  A  contains  the  same  data;  however,  the  data  has  been  formed 
into  dimensionless  groups  for  the  purpose  of  analyzing  and  comparing 
resistance  obtained  in  this  series  with  earlier  series  results  and 
eliminating  the  effects  of  variations  in  thickness  from  test  to  test. 

The  third  section  of  Appendix  A  contains  the  data  scaled-up  using  the 
scaling  laws  in  Table  2.2  to  reach  values  of  resistance  and  velocity, 
etc., for  each  individual  data  point  as  they  would  pertain  to  the 
actual  tug.  Table  3.1  contains  conversion  factors  from  English  to  SI  Units. 

The  model  data  is  listed  in  SI  units.  This  system  was  commonly 
known  as  the  metre  kilogram  second  system  up  to  1960.  The  resistance  is 
thus  measured  in  newtons,  the  velocity  in  millimetres  per  second,  the 
thickness  of  the  ice  in  millimetres,  the  strength  of  the  ice  in  kilo- 
pascals,  the  pressure  in  the  manifolds  is  in  pascals, and  the  temperature 
measured  in  the  manifolds  in  degrees  Celsius.  The  volume  flow  rate  is 
given  in  cubic  metres  per  hour,  while  the  power  is  given  in  watts.  The 
volume  flow  has  been  calculated  for  standard  conditions  using  Equation 
2.1,  which  is  repeated  here. 
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The  compressor  power  necessary  to  obtain  this  volume  flow  rate  is  cal¬ 
culated  using  a  compressor  efficiency  of  0.7  and  the  equation  for 
isentropic  compression  of  air  is  as  follows: 


kp  xVi 

'?2 

k- 1 

w  -  k=r[ 

1  Pi 

R  -  1 

(3.1) 


where : 


y  —  work 

k  =  ratio  of  specific  heats,  1.4  for  air 
Pj  =  initial  pressure 
p2  =  outlet  pressure 
vx  -  inlet  volume 
R  =  gas  constant  for  air 


Pressure  at  the  air  compressor  on  the  tug  itself  has  been  calculated 
from  the  volume  flow  rate  using  the  actual  piping  system  performance 
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TABLE  3.1  CONVERSION  FACTORS 


To  Convert  From 

Jo 

Multiply  By 

long  tons 

kN 

9.962 

ft 

m 

0.3048 

knots 

m/s 

0.5144 

lbf/in2 

kPa 

6.895 

ft3/min 

m3/s 

4.719xl0'4 

HP 

watts 

746.0 

1 I 


as  calculated  in  the  design  manual  which  is  Reference  2.  This  manual 
contains  performance  curves  of  volume  flow  rate  versus  pressure  at 
the  compressor. 

The  dimensionless  groups  which  were  used  in  this  program  were 
those  which  have  been  developed  using  long  standing  practice  at 
ARCTEC.  The  quantity  of  data  available  at  any  single  condition  of  the 
hull  or  bubbler  system  is  too  small  to  allow  the  empirical  formation  of 
different  dimensionless  groups.  One  main  purpose  of  dimensionless 
groups  is  to  reduce  the  number  of  variables  present  such  that  the 
analysis  can  be  more  efficiently  conducted.  The  single  variable  which 
effects  the  performance  most  is  the  thickness  of  the  ice;  therefore, 
each  dimensionless  group  contains  the  thickness  of  the  ice  raised  to 
some  power  within  it.  Figure  3.1  shows  a  plot  of  the  dimensionless 
resistance  versus  the  product  of  dimensionless  strength  and  dimension¬ 
less  velocity.  These  two  terms  were  found  to  be  the  most  important 
in  explaining  the  variation  in  resistance  found  in  the  earlier  WYTM 
tests.  All  of  the  data  obtained  during  this  test  program  in  solid 
ice  is  plotted  in  Figure  3.1  using  these  two  dimensionless  variables. 

A  look  at  the  data  indicates  that  there  is  a  large  amount  of  scatter 
present  and  that  no  direct  conclusions  are  possible  by  simply  looking 
at  the  way  the  data  falls.  In  brash  ice,  previous  experience  leads 
to  the  formation  of  two  different  dimensionless  groups  to  explain 
the  variation  in  resistance  which  is  seen.  The  first  dimensionless 
group  or  dimensionless  brash  resistance  is  the  resistance  measured 
divided  by  the  weight  density  of  ice  times  the  beam  of  the  model  and 
thickness  of  the  brash  field  squared.  The  dimensionless  velocity 
squared  is  the  term  against  which  this  resistance  is  plotted.  It  can 
also  be  seen  from  an  examination  of  the  plot  in  Figure  3.2  that  a 
readily  made  conclusion  is  not  possible  on  the  performance  of  the  tug 
in  brash  due  to  the  large  scatter  in  the  data. 

We  can,  however,  compare  the  results  which  were  obtained  in 
this  series  of  tests  without  the  bubbler  system  operating  with  those 
which  were  obtained  using  the  sanded  hull  form  during  the  earlier 
tests  of  the  140-FT.  WYTM.  Data  which  falls  within  the  range  of  data 
gathered  during  this  series  is  plotted  on  Figure  3.3.  A  separate 
regression  equation  was  developed  for  the  tests  conducted  earlier  in 
only  1-1/2  feet  of  ice.  This  line  is  also  plotted  on  Figure  3.3. 

Along  with  the  data  collected  during  these  tests,  we  have  plotted  a 
regression  of  the  performance  of  the  tug  without  the  bubbler  system 
also.  From  this  it  can  be  seen  that  in  spite  of  the  fact  that  a 
lower  hull-ice  friction  was  measured  in  this  series,  the  resistance 
was  somewhat  higher  at  all  speeds  tested  than  for  the  earlier  tests. 
However,  we  can  also  see  that  the  scatter  of  the  data  is  such  that 
all  of  the  data  points  collected  both  before  this  series  and  during 
this  series  fall  within  the  same  band,  and  thus  represent  the  same 
group.  The  error  or  mismeasurement  of  the  resistance  appears  to  vary 
from  approximately  10%  at  zero  speed  to  33%  at  5  knots. 
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3.2  Friction  Data 


The  hull-ice  friction  was  measured  at  four  locations  on  the 
hull  of  the  140-FT.  WYTM  Model.  The  friction  was  measured  many 
times  at  each  of  the  four  locations  and  this  data  is  recorded  in 
Appendix  B.  The  average  of  all  the  friction  tests  was  0.232. 

3.3  Ice  Data 

The  properties  of  the  model  ice  are  contained  in  Appendix  C. 
The  strength  given  in  this  particular  tabulation  is  the  average  of 
the  strength  for  all  data  points  obtained  in  solid  ice  during  a  given 
day,  thus,  it  may  not  correspond  exactly  to  the  strength  recorded  for 
individual  data  points.  In  the  analysis,  individual  strengths  were 
used  for  each  data  point.  In  the  brash  data,  it  should  be  noted  that 
the  density  of  ice  to  water  within  the  brash  sheet  was  approximately 
0.65. 


4.  ANALYSIS 


4.1  Tug  Performance  in  Open  Water 

Thrust  curves  for  the  actual  propeller  to  be  utilized  and 
the  actual  hull  form  for  the  WYTM  were  obtained  using  the  data  con¬ 
tained  in  Reference  5.  The  1-t  value  of  0.81  and  the  l-lfy  value 
of  0.745  tested  at  5  knots  were  assumed  to  hold  constant  down  to 
the  bollard.  The  propeller  thrust  was  calculated  from  the  Equation: 

K  pJl2Dk0-t) 

Tp  =  ~ — -  =  3,66  kt  in  1ong  tons  (4J) 

While  the  ship  velocity  at  various  J-factors  was  calculated  using 
this  following  equation: 


vk  ~  oI'8"(T-¥y]'  =  5-76nJ  in  knots  (4-2) 


Open  water  resistance  was  calculated  using  the  following  equations 
and  the  values  for  estimated  horsepower  tabulated  in  the  David  Taylor 
Model  Basin  tests. 
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The  curves  were  constructed  for  five  values  of  constant  shaft  horse¬ 
power  varying  from  the  bollard  to  about  5  knots.  The  results  of 
these  calculations  are  presented  in  Figure  4.1  which  is  a  re¬ 
presentation  of  the  actual  performance  of  the  WYTM  with  its  installed 
power  for  intermediate  values  of  horsepower  up  to  full  power  on  the 
shaft.  Also,  plotted  on  this  figure  are  the  open  water  resistance  and 
three  additional  power  levels  which  are  sketched  in  above  the  capa¬ 
bility  of  the  installed  propulsion  plant,  for  the  purpose  of  evaluating 
some  of  the  solid  ice  resistance  data  at  5  knots  which  was  much  higher 
than  the  capabilities  of  the  present  ship. 


4.2  Continuous  Motion  Through  Solid  Ice 

In  this  section  we  will  discuss  the  results  which  were  ob¬ 
tained  for  towing  the  ship  at  constant  speed  through  solid  ice  with  and 
without  the  air  bubbler  system  in  operation.  The  data  from  Appendix  A 
was  placed  into  a  file  on  the  computer  and  submitted  to  a  forced  fit  to 
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the  following  two  equations: 
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The  results  of  this  regression  analysis  are  tabulated  in  Table  4.1 
which  contains  the  regression  coefficients  for  resistance  in  solid 
ice.  The  first  column  in  this  table  contains  the  setting  names 
which  pertain  to  the  settings  or  division  of  air  flow  to  the  bow  and 
the  midships  manifold  given  in  Table  2.3.  The  second  column  is  the 
volume  flow  in  standard  cubic  feet  per  minute,  full-scale  corrected 
for  the  temperature  and  pressure  of  the  gas  utilized  during  the 
tests.  The  third  column  is  the  full-scale  compressor  power  required 
based  on  the  existing  WYTM  piping  system  and  the  pressure  drop 
associated  with  that  system  as  calculated  in  the  design  report. 
Reference  2.  The  fourth,  fifth,  and  sixth  columns  in  this  table  are 
the  mean  values  of  the  variables  which  were  subjected  to  regression 
for  each  of  the  conditions  tested.  The  regression  solution  always 
forces  the  regression  line  through  the  mean  value  of  all  the  data, 
therefore,  the  line  Equation  4.4  must  pass  through  the  mean  value  of 
R/p ugh2  at  the  mean  value  of: 


Jgh) 

while  equation  4.5  must  utilize  the  mean  value  of  R/pwgBh2  and  the 
mean  value  of  v/Jgh.  Computer  print-outs  for  these  fits  are  found 
in  Appendix  D.  The  coefficients  that  were  determined  by  the  regression 
solution  are  then  tabulated  in  the  next  four  columns  listed  CQ ,  C\3 
Cz  and  C3.  The  number  of  data  points  which  were  available  for  each 
regression  are  tabulated  in  the  last  column.  Because  of  the  small 
number  of  data  points,  the  value  for  the  multiple  correlation  coefficient 
for  these  regression  equations  was  quite  poor.  However,  the  fit  to 
the  first  equation,  Equation  4.4,  was  much  better  than  the  fit  to 
Equation  4.5.  Also,  as  can  be  seen  by  a  cursory  examination  of  the 
coefficients,  the  agreement  between  the  different  flow  rates  is  much 
better  when  the  strength  effect  is  taken  into  consideration,  as  it  is 
in  Figure  4. 4, than  when  it  is  not  as  in  Equation  4.5. 
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The  coefficients  C0  and  c,  were  then  plotted  against  volume 
flow  on  Figure  4.2.  The  slopes  for  Q3a,  Q2,  and  Q3  are  omitted  from 
the  figure  because  they  appeared  to  be  abnormally  high.  A  line  passed 
through  the  slopes  for  Q0  and  Q1  (those  points  containing  the  most 
data)  passes  near  the  slopes  of  Qh„  and  Qh.  This  line  was  used  to 
obtain  "faired"  values  of  slope,  C\,  for  the  remainder  of  the  analysis. 
The  corresponding  value  of  intercept,  c0,  to  pass  a  line  through  the 
mean  of  each  flow  setting  data  set  is  then  calculated  in  the  next 
column  of  Table  4.2.  The  second  intercept  values  are  also  plotted 
in  the  lower  half  of  Figure  4.2  and  a  line  drawn  through  them.  The 
values  along  the  line  were  then  read  off  and  a  smooth  family  of  re¬ 
sistance  curves  is  thus  obtained  such  that  the  mean  for  the  family 
goes  through  the  mean  of  the  raw  data.  This  family  of  dimensionless 
resistance  curves  is  plotted  in  Figure  4.3.  Also  plotted  in  Figure 
4.3  are  the  thrust  curves  for  the  tug  when  the  resistance  curves 
represent  1-1/2  feet  of  126.5  psi  ice.  The  dimensionless  coefficients 
could  be  used  to  calculate  the  resistance  in  any  other 
thickness  of  ice  desired.  The  power  required  for  the  tug  to  proceeu 
through  1-1/2  feet  of  ice  was  then  evaluated  at  1 ,  3  and  5  knots.  The 
results  of  this  evaluation  are  contained  in  Figures  4.4  through  4.6. 
This  is  plotted  as  the  required  horsepower  at  this  speed  through 
1-1/2  feet  of  126.5  psi  ice.  The  value  of  126.5  psi  is  the  mean  value 
of  ice  strength  encountered  during  the  tests  of  the  icebreaker 
MACKINAW  on  the  Great  Lakes.  The  solid  line  on  each  of  these  curves 
is  the  compressor  horsepower  required  to  produce  the  volume  flow 
through  the  existing  piping  system  on  the  WYTM.  If  the  piping  system 
were  resized  and  optimized  for  each  particular  flow  rate,  the  pressure 
drop  would  remain  constant  and  this  line  would  become  straight  as  shown 
by  the  dashed  line  on  Figure  4.4.  Since  7,500  scfm  is  the  design 
point,  a  straight  line  through  the  7,500  scfm  power  point  and  the 
origin  would  be  the  characteristic  of  the  compressor  power  curve  under 
these  conditions.  The  shaft  horsepower  curve  is  extracted  from  the 
predictor  equations  and  their  intercept  with  the  propeller  character¬ 
istics.  The  total  power  is  then  the  sum  of  compressor  power  and  shaft 
horsepower.  The  lowest  point  on  the  total  power  curve  corresponds  to 
the  point  of  optimum  operation  of  the  bubbler  system  for  that  partic¬ 
ular  speed. 


4.3  Continuous  Motion  Through  Brash  Ice 

The  continuous  motion  through  brash  ice  data  contained  in 
Appendix  A  was  formed  into  dimensionless  groups  corresponding  to 
the  groups  in  the  following  equation: 
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TABLE  4.2  PREDICTOR  FOR  SOLID  ICE  \ 


Setting 

Name 

Vol . 
ft3/m 

^estimated 

Co* 

^estimated 

Qo 

0 

.0307 

7.82 

7.82 

1,000 

.0314 

8.80 

7.63 

2,000 

.0321 

7.68 

7.45 

Qi 

4,000 

.0334 

7.14 

7.19 

Qz 

8,000 

.0361 

6.63 

6.50 

Q  3 

12,250 

.0390 

5.62 

5.80 

Q4 

16,800 

.0420 

4.11 

5.05 
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Where: 


R  =  resistance  to  movement 

p.  =  mass  density  of  ice 

g  =  gravitational  constant 

B  -  beam  at  the  waterline 
h  =  ice  or  brash  thickness 
v  =  velocity  of  the  ship 


The  data  for  each  separate  flow  rate  was  submitted  to  a  forced  re¬ 
gression  fit  and  the  results  of  this  are  tabulated  in  Table  4.3  As 
in  Table  4.2,  the  first  five  columns  correspond  to  the  settings  em¬ 
ployed  during  the  test  program,  the  volume  flow  rate  in  standard 
cubic  feet  per  minute  at  full-scale  for  these  conditions,  the  com¬ 
pressor  power  at  full-scale  for  these  conditions,  the  mean  value  of 
the  dimensionless  brash  resistance  and  the  mean  value  of  the  square 
of  dimensionless  velocity.  The  two  coefficients  obtained  through 
regression  analysis  of  this  data  follow  with  the  last  column  con¬ 
taining  the  number  of  points  used  to  determine  these  coefficients. 

The  same  general  process  employed  in  analyzing  the  continuous  motion 
through  solid  ice  was  again  used  to  determine  predictor  for  contin¬ 
uous  motion  through  brash  ice.  First,  a  plot  was  made  of  the  slopes 
for  each  of  the  equations  as  computed.  It  is  very  difficult  to  see 
where  a  line  should  be  drawn  through  these  points,  therefore,  the 
first  three  points  and  last  four  points  are  each  averaged  and  plotted 
as  triangles  on  Figure  4.7.  A  line  drawn  through  these  two  mean 
points  is  used  to  develop  a  family  of  curves  having  changing  slope 
and  also  fitting  the  mean  of  all  the  data.  These  new  slopes  were 
employed  as  shown  in  Table  4.4  to  calculate  a  new  intercept  value  for 
these  conditions.  The  C4  and  Cs  estimated  values  represent  a  smooth 
family  of  curves  which  can  be  used  to  predict  the  performance  of  the 
140-FT.  WYTM  through  brash.  One  serious  problem  with  the  data  from 
the  brash  tests,  is  that  the  data  was  obtained  in  3  ft  of  ice,  well 
below  the  maximum  capabilities  of  the  WYTM.  These  predictor  equations 
from  Table  4.4  were  employed  to  develop  the  resistance  in  6  ft  of 
brash  from  data  obtained  in  3  ft  of  brash.  This  has  its  limitations 
in  that  it  must  be  assumed  that  the  extrapolation  is  correct  from  3  ft 
to  6  ft.  On  the  other  hand,  it  gives  us  a  better  idea  of  the  perfor¬ 
mance  of  the  tug  nearer  its  limiting  conditions.  The  most  likely 
disadvantage  of  this  method  is  that  it  may  emphasize  the  effect  of 
the  air  bubbler  system  as  the  dimensionless  equation  is  extrapolated 
to  twice  the  test  thickness.  The  predictor  equations  are  plotted  on 
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TABLE  4.3  REGRESSION  OF  BRASH  RESISTANCE  DATA 


Compressor 

Name  Vol .  Power 
ft3/m  HP 


R 


p^Bh2- 

mean 


(£)' 

mean 


C  4 


#  Of 
Points 


TABLE  4.4  BRASH  PREDICTOR  EQUATION 


Name 

Vol . 

^estimated 

<?4* 

^“estimated 

«0 

0 

0.100 

0.524 

0.600 

B 

1,000 

0.1 30 1 

0.570 

0.582 

«.  44 

2,000 

0.162 

0.534 

0.562 

4,000 

0.222 

0.633 

0.522 

«2 

8,000 

0.345 

0.535 

0.444 

10,125 

0 . 41 3  2 

«3 

12,250 

0.476 

0.353 

0.362 

16,800 

0.616 

0.279 

0.275 

1  (Cs)0  +  (C s)0  +  (^s) 


34 


3 


i(<7s)«,  +  (Cs)«2  +  (Cs)e,  +  (Ci)e. 


4 
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both  the  dimensionless  basis  and  the  six  feet  of  brash  ice  in  Figure 
4.8.  The  ship  performance  curves,  of  course,  only  apply  to  the 
dimensional  axes.  Intersections  of  the  resistance  curves  and  the 
thrust  curves  are  operating  points  at  which  the  tug  will  proceed  in 
a  continuous  fashion  through  6  ft  of  brash.  This  figure  was  used 
to  evaluate  the  performance  of  the  bubbler  system  at  1 ,  3  and  5  knots 
through  6  ft  of  brash  ice.  This  information  is  then  plotted  in  Figures 
4.9  through  4.11.  The  plots  were  made  in  the  same  fashion  as  those 
for  the  uniform  sheet  ice  and  again  the  optimum  operating  condition 
is  the  minimum  total  power  point  as  shown  on  these  curves.  One  thing 
that  we  do  see  in  the  brash  ice  is  that  there  are  occasions  when 
more  air  flow  is  required  than  the  presently  installed  compressor  can 
provide. 


4.4  Starting  Performance 


In  addition  to  the  towed  resistance  tests,  a  series  of 
starting  resistance  tests  were  performed  in  both  solid  and  brash 
ice.  This  was  the  first  time  that  tests  of  this  sort  have  been 
performed  at  ARCTEC  and  a  number  of  initial  problems  were  encountered. 
The  time  required  to  get  the  towing  apparatus  started  allowed  creep 
to  take  place  in  the  ice  sheet  and  very  low  values  of  resistance  were 
obtained  during  the  first  two  days  of  the  tests.  On  the  last  two 
days,  however,  the  procedure  was  refined  such  that  the  ship  was 
forced  by  hand  into  the  ice  sheet  and  quickly  attached  to  the  towing 
apparatus  for  the  very  slow  speed  starting  force  determination. 

Very  consistent  results  were  obtained.  The  best  of  these  results 
are  plotted  in  Figure  4.12.  Starting  resistance  is  plotted  as  a 
function  of  the  bubbler  flow  for  three  particular  ice  conditions  in 
which  tests  were  actually  conducted.  Also  plotted  on  this  figure, 
as  dotted  lines, are  the  propeller  thrust  of  the  installed  propulsion 
system.  Using  the  propeller  thrust  curves,  the  necessary  shaft 
horsepower  to  start  the  tug  moving  was  then  calculated  and  plotted 
in  Figure  4.13.  When  this  shaft  horsepower  is  added  to  the  compressor 
power,  a  total  power  curve  is  produced  which  is  similar  to  those  pro¬ 
duced  for  continuous  motion  in  solid  ice  and  brash  earlier.  For  the 
starting  resistance,  it  appears  that  a  much  lower  bubbler  flow  is 
necessary  to  achieve  optimum  results.  A  flow  of  approximately  4,000 
scfm  in  3.6  feet  of  brash  gives  the  best  performance.  Only  slightly 
higher  flows  are  required  for  the  two  sheet  ice  tests  which  were  per¬ 
formed.  There  should  be  caution  applied  to  these  performance  curves 
in  that  the  ice  was  a  great  deal  weaker  than  the  target  value  of 
126.5  psi.  The  lower  curve  in  fact,  corresponds  more  to  melting  sea 
ice  whereas  the  upper  curve  is  intermediate  betweln  sea  ice  and  fresh 


£)i m  3ns  ion/bsj  ^ 


igura4-Q:Tug  */ith  &ubb/»r  in  6-0  ft  &mstt 


/to/~j»joaL/er-  /r  ed  *  to 


Figur*  4.9: &**bb/»r  PerForntanee  at  f  Knot  in  6.0  Ft  &ras A 


ifciaiki  >ih  >*  >i«ifrgai» . 


4-18 


Morse  fio»rer  Megu/reM  <  fO 


Figure  4.10  •  Bubbler  Performance  at  3  Knots  in  6.0  ft  &rash 
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Figure  4.  H  (  Bubbler  Performance  a*  S  Knots  in  6.0  ft  Brash 
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water  ice.  The  data  points  from  these  three  curves  is  plotted 
on  a  dimensionless  basis  in  Figure  4.14.  This  shows  the  actual 
percentage  reduction  in  starting  resistance  due  to  the  installation 
of  an  air  bubbler  system.  This  also  clearly  shows  that  the  design 
flow  is  sufficient  to  produce  nearly  minimal  starting  resistance. 


4.5  Analysis  of  the  Photographic  Data 

The  underwater  movies  which  were  taken  clearly  show  the 
change  in  performance  of  the  air  bubbler  system  with  speed  and  flow 
rate.  The  nozzles  on  the  bow  manifold  do  not  discharge  downward  so 
that  seachests  located  aft  of  the  bow  manifold  would  not  appear  to 
be  subject  to  injestion  of  air.  The  movies  also  clearly  show  that 
the  train  of  bubbles  from  the  nozzles  rising  at  a  constant  speed  is 
swept  further  and  further  aft  as  the  speed  of  the  tug  increases. 
Therefore,  the  section  of  the  bow  subjected  to  the  heaviest  amount 
of  breaking  is  never  covered  with  air  and  the  bow  shoulders  are 
covered  with  air  less  as  speed  increases.  There  is  a  great  deal  more 
film  available  for  the  condition  of  operation  in  sheet  ice  than  there 
is  available  to  study  operation  of  the  bubbler  system  in  brash.  In 
sheet  ice  the  bubbles  are  trapped  between  the  surface  of  the  ice  and 
the  hull  and  seem  to  effectively  keep  the  ice  away  from  the  hull.  In 
brash,  because  of  the  broken  coating,  there  appears  to  be  some  leak¬ 
age  of  the  air  out  from  the  side  of  the  ship  and  the  uneven  surface 
requires  more  air.  If  possible,  the  technical  film  to  examine  how 
the  bubbles  travel  along  the  hull  should  be  viewed. 


3u  hbfar- 
ur/thouf  3 ubbJ 


□i 


Figure  4.  f4  1  Qiman-a  ion/eaa  'Starting  /?es /a+artcm. 
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5.  CONCLUSIONS  AND  RECOMMENDATIONS 


5.1  Conclusions 


The  bubbler  system  as  designed  at  7,500  scfm  is  effective 
in  reducing  the  resistance  of  the  icebreaking  tug  and 
also  in  reducing  its  fuel  consumption.  As  can  be  seen 
in  the  performance  curves  plotted  in  section  4,  in  every 
case  the  total  horsepower  required  to  proceed  through 
a  given  ice  condition  was  reduced  due  to  use  of  the  air 
bubbler  system.  The  system  is  most  effective  at  low 
velocities  and  in  brash  ice.  For  starting  resistance 
in  brash  ice, the  resistance  will  only  be  60%  of  that 
without  a  bubbler  system  when  the  bubbler  system  is 
operated  at  its  designed  conditions.  For  starting 
conditions  in  solid  ice, the  resistance  is  reduced  by 
20%  with  operation  of  a  bubbler  system  at  its  designed 
conditions. 

A  good  distribution  of  the  air  coating  is  obtained  along 
the  shoulders  and  side  of  the  tug  by  the  air  bubble  train. 
There  appears  to  be  no  need  to  place  intermediate  nozzles 
in  order  to  get  a  satisfactory  distribution  of  air. 


Most  of  the  bow  is  not  covered  by  air  bubbles.  In 
future  designs  more  effort  needs  to  be  given  to  the 
possibility  of  locating  air  nozzles  further  forward. 

As  was  shown  in  the  coating  tests  completed  prior  to 
this  program,  the  area  immediately  adjacent  to  the 
stem  is  that  contributing  the  most  to  icebreaking  re¬ 
sistance.  It  therefore  seems  logical  that  to  be  most 
effective  an  air  bubbler  system  should  coat  this  par¬ 
ticular  area. 

5.2  Recorrcnendations 

As  soon  as  the  full-scale  WYTM  becomes  available  in  its 
initial  encounter  with  ice,  a  series  of  full-scale  tests 
should  be  conducted  to  obtain  at  least  a  few  data  points 
which  could  be  used  to  compare  with  these  model  scale 
results.  This  would  then  improve  the  confidence  in  the 
effects  of  various  flow  rates  on  performance  as  was 
determined  in  this  series. 
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2.  The  use  of  additional  nozzles  forward  of  those  used  in 
this  program  should  be  examined,  particularly  in 
future  designs.  It  would  be  possible  to  equip  the  model 
of  the  140  FT.WYTM  tug  with  additional  nozzles  further 
forward  to  study  this  particular  possibility  before 

any  further  design  work  is  accomplished. 

3.  A  study  should  be  conducted  of  the  trade-off  between 
employment  of  a  air  bubbler  system  and  application  of 
a  low  friction  underwater  coating.  It  is  known  now 
that  both  of  these  systems  will  improve  the  performance 
of  icebreaking  ships  but  the  interaction  between  the 
two  is  still  not  known.  For  example, we  do  not  know  if 
the  bubbler  system  mechanism  for  reducing  icebreaking 
resistance  is  the  same  as  that  of  the  low  friction 
coating  :in  that  forming  a  film  of  air  and  water  between 
the  ice  and  the  hull  it  reduces  hull/ice  friction. 

This  question  might  be  answered  by  coating  the 

140  FT.  WYTM  with  a  smooth  paint  surface  and  repeating 
a  few  tests  with  the  bubbler  to  see  if  the  benefits 
are  still  obtained  with  the  smooth  hull. 


5-2 


6.  REFERENCES 


1.  Lecourt,  E.J.,  "Icebreaking  Model  Tests  of  the  140-F00T  WYTM," 
ARCTEC  Report  202C-2,  Contract  D0T-CG-50383-A,  September  1975. 

2.  Lecourt,  E.J.,  "Air  Bubbler  System  Design  for  the  140-F00T  WYTM," 
ARCTEC  Report  202C-3,  Contract  D0T-CG-50383-A,  January  1976. 

3.  ARCTEC,  Incorporated  "Test  Plan  for  Model  Tests  in  Ice  to  Confirm 
Effectiveness  of  the  140-F00T  WYTM  Air  Bubbler  System,"  ARCTEC 
Report  354C-1 ,  Contract  D0T-CG-64.242-A,  December  1976. 

4.  ARCTEC,  Incorporated  "Proposal -Study  to  Conduct  Model  Testing 
in  Ice  to  Confirm  Effectiveness  of  140-F00T  WYTM  Air  Bubbler 
System,"  Solicitation  CG-64-242A,  August  1976. 

5.  Remmers,  K.D.  and  R.  Hecker,  "Powering  Predictions  for  the  United 
States  Coast  Guard  140-F00T  WYTM  [Icebreaking  Tugboat]  (Model 
5336-Propeller  4657)"  NSRDC  Report  SPD-223-18,  July  1975. 

6.  Baumeister,  Theodore,  editor,  Mark's  Mechanical  Engineer's  Hand¬ 
book,  Sixth  Edition,  McGraw  Hill,  New  York,  1958,  Chapter  14. 

7.  Lewis,  J.W.,  "Ship  Model  Ice  Resistance  Experiments,"  Technical 
'  Report  0029,  Contract  No.  D0T-CG-13189-A,  ARCTEC,  Incorporated, 

October,  1972. 

8.  Kobus,  Helmut  E.,  "Analysis  of  the  Flow  Induced  by  Air  Bubbler 
Systems,"  ASCE,  Proceedings  of  the  Eleventh  Conference  on  Coastal 
Engineering,  London,  England,  September,  1968. 


APPENDIX  A 


RESISTANCE  DATA 


TABLE  A.l  Model  Data 

Solid  Ice 

Brash 

Starting 


TABLE  A. 2  Dimensionless  Data 

Solid  Ice 
Brash 


TABLE  A. 3  Data  in  Full-Scale  Units 
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TABLE  A-l .  MODEL  DATA 
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+  Off  Scale,  estimated  18.7  kPa. 


TABLE  A-l.  (Continued) 
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TABLE  A-l.  (Continued) 
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Deleted  from  regression  analysis  due  to  wide  fluctuation  in  load. 


_Eh  -•-» 

Hi  « 


00  00  *3-  CO 

VOUOr—  O  I —  > —  00  CO  i—  00  CO 

i—  00  CNJ  i—  r-CVJCV  MOWOMM 

icjcvjooni  i  cuon  i  upincoiCN  i  i  i  i 

I  ••••II  ...I . I  I  I  I 

n  r-f-  I-  i —  i —  CO 


Eh  ^ 
O  m 
E 


CD  o  CTt  r—  CD  O  CM  CD  •—  CD  O  CM 


co  »—  co  i 

•  •  •  •  I 

O  CM  CM 


f—  O  O  O  I 
i  I 


r-  o  O  O  I 

I  I 


I—  O  O  O  I 
I  I 


CM  O  O  O  I 


CO  CM  CM 

r-v  cmn  j 

•  +  •  •  I 

if)  lo  rs  ^ 

<3- 

cm  cm  vo  cd 

r—  CO  I 

•  -  •  I 

00  CM  f— 


r-s.  r—  ^ 

o>  o  i 

.  •  •  •  i 

o  r— 


oo  m  o\  m 

^  O)  N  I 
«  «  •  «  | 
LD  CM  •— 


I  1 

I  I  CO- 
CM 


VO  »—  ^  I  CM  lO  lO  »—  ON  I 

•  •  . . I 

r—  CM  »“  CM  r— 


10  0  0  0=0  I 

I  I 


I  o  o  o  o  -  O  I 
I  I 


10  0  0  0=0  I 


10  0  00=0  I 
I  I 


«—  ^  00 
CM  CM  *3“ 


r—  CO  00 

cm  rr 


VO  cm  in 

CM  i—  CM 


^ 

CM  rs  CM 
*  ♦  ♦ 
CM  r—  CM 


in  vo  o>  in 

i  <3- cd  cd  i  i  i  i 

i  •  •  ■  •  =  *iiii 

CM  CM  CD  r-  CM 


U)(SJOO)  'Tf 
I  CMiOO)^-  .  CO  I 

I  •  •  •  •  =  •  I 

r-  CM  r-  00 


CM  CM  CM  CM 

I  VO  r—  VO  _  I —  I 

I  .  #  •  .  -  •  I 

r-  r—  CM  CM  CM 


cd  cm  cm  a>  cr> 

*  .  CD  I 

I  •  *  •  *  -  •  I 

,-r-CMr-  CM 


t  ;  i  i 
i  •  i  i 


o>o»wcnounoio«t^-oi^-^'0>u:oo«3^,fl-'t 

>—  i—  r-  f—  r—  r-<T>i—  OOr—  00<—  (^-CDCTlOOOO 


VO 

o\ 

CD 

in 

• 

_  *  - 

CM 

CM 

-  VM 

CM 

CM 

vo 

U5  § 

o  4- 

55  4>  ** 
5  u  •- 
r~  *->  -a 


ouflr-inmoTj-ai^ai^t'CMOooo^’OOOivrs 

icHONNr-tncncvjmioi-coivacoM'H'-O't'- 

CMN^'toooowcicoMNcvjomn^'CO^rsm 

-C  T3  T3  T3  -C 

(/I  i—  *r-  (/) 

IQ  ^  ^  r—  ^  »— - ^ 

1-0  0  O  1- 

co  in  ts>  to  co 


05  ^  05  05  "3; 

ocHJ’^-^-oo^cnrHom^F-jrncHoooo 

cycycycyQrQpQfOocyQpcycafQrcycycycyQi’cyQf 


■NPm’lfMOrsWWOr-Cvm^'lflOPsCOOOr- 

r—  *—  r—  i—  •—  i—  r—  r—  f—  I—  C\J  C\J 


\  i 

i  j 


TABLE  A-2  DIMENSIONLESS  DATA 


SHEET  ICE 


# 

Q 

R 

_v _ 

a 

a'  x  v' 

WBh' 

Sgh 

1 

5o 

12.4 

.206 

297 

61.3 

2 

16.0 

.644 

291 

187 

3 

1 

! 

19.6 

1.167 

308 

360 

4 

Qtt 

6.78 

.167 

291 

48.7 

5 

13.1 

.639 

286 

183 

6 

1 

16.3 

1.132 

299 

338 

7 

10.7 

.713 

175 

125 

8 

1 

16.1 

1.192 

185 

221 

9 

£ 

5.63 

.172 

165 

28.3 

10 

17.2 

1.156 

182 

211 

11 

< 

h 

6.40 

.172 

175 

30.0 

12 

1 

10.6 

.645 

173 

111 

13 

8.83 

.174 

245 

42.7 

14 

10.51 

.602 

237 

142.6 

15 

13.66 

1.032 

228 

235.2 

16 

£ 

o 

8.54 

.162 

240 

38.8 

17 

9.44 

.692 

248 

171.8 

18 

13.33 

1.068 

233 

249.1 

19 

«3 

8.31 

.727 

148 

107.5 

20 

1 

12.85 

1.077 

143 

154.4 

21 

4.10 

.164 

132 

21.6 

22 

12.17 

1.151 

146 

167.5 

23 

4.46 

.169 

143 

24.2 

24 

i 

8.59 

.632 

137 

86.8 

25 

10.21 

.642 

84 

53.7 

26 

12.31 

1.009 

83 

83.6 

27 

7.49 

.160 

83 

13.3 

28 

i 

11.18 

1.018 

84 

85.9 

29 

£ 

u 

4.28 

.154 

83 

12.8 

30 

1 

8.09 

.575 

82 

47.4 

31 

5.63 

.170 

70 

11.8 

32 

il 

7.95 

.619 

72 

44.5 

33 

9.87 

.998 

67 

67.1 

34 

4. 

5.18 

.162 

70 

11.3 

35 

9.01 

.581 

70 

40.6 

36 

1 

10.44 

1.013 

69 

70.2 

37 

fl2 

14.3 

.655 

157 

103 

38 

i 

15.4 

1.189 

148 

176 

39 

10.0 

.166 

144 

23.8 

40 

14.7 

1.056 

147 

155 

A-8 


TABLE  A-2  DIMENSIONLESS  DATA  (Continued) 

SHEET  ICE 


# 

Q 

R 

V 

Cf 

<J'  X  V 

WBhi 

Jgh 

41 

Q: 

6.5 

.162 

144 

23.2 

42 

i 

12.0 

.627 

144 

90.0 

43 

Qhb 

9.3 

.176 

123 

21.7 

44 

11.1 

.610 

127 

77.8 

45 

l 

12.8 

1.151 

117 

135 

46 

8.5 

.172 

118 

20.2 

47 

1 

10.9 

.585 

117 

68.5 

48 

1 

13.6 

1.032 

119 

122 

p„  •  1005  Kg/m3 
g  =  9.80  m/sec? 
3  =  0.434  m. 
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TABLE  A- 2  DIMENSIONLESS  DATA  (Continued) 


BRASH  ICE 


# 

Q 

R 

V 

R 

V2 

gh 

P^h1 

Jgh 

pTpF" 

1 

Qi 

.423 

.455 

.450 

.207 

2 

.844 

.803 

.898 

.645 

3 

Q  3 

.267 

.125 

.284 

.016 

4 

.465 

.578 

.495 

.334 

5 

.480 

.108 

.511 

.012 

6 

J 

.421 

.410 

.448 

.168 

7 

Q  i 

.501 

.117 

.533 

.014 

8 

.637 

.428 

.678 

.183 

9 

.981 

.721 

1.044 

.520 

10 

Q  o 

.530 

.171 

.564 

.029 

11 

.383 

.458 

.407 

.210 

12 

.553 

.647 

.588 

.419 

13 

«3 

.421 

.426 

.448 

.181 

14 

J 

.671 

.748 

.714 

.559 

15 

.193 

.123 

.205 

.015 

16 

.325 

.677 

.346 

.458 

17 

«a 

.324 

.150 

.345 

.023 

18 

.459 

.493 

.488 

.243 

19 

.602 

.446 

.640 

.199 

20 

1 

! 

.844 

.826 

.898 

.682 

21 

.722 

.161 

.768 

.026 

22 

1 

.899 

.729 

.956 

.531 

23 

h 

.903 

.129 

.961 

.017 

24 

J 

.587 

.459 

.624 

.211 

25 

•u 

.395 

.126 

.420 

.016 

26 

.438 

.450 

.466 

.202 

27 

.510 

.683 

.542 

.466 

28 

.592 

.119 

.630 

.014 

29 

.582 

.405 

.619 

.164 

30 

.689 

.740 

.733 

.548 

31 

«a 

.43/ 

.466 

.465 

.217 

32 

* 

.738 

.769 

.785 

.591 

33 

934 

.773 

.264 

.822 

.070 

34 

J 

.677 

.782 

.720 

.611 

35 

1 1 

.533 

.248 

.567 

.061 

36 

. 

.312 

.448 

.332 

.201 

37 

< 

■43 

.656 

.282 

.698 

.079 

38 

.547 

.453 

.582 

.205 

39 

.614 

.763 

.653 

.582 

40 

«o 

.560 

.237 

.596 

.056 

A-10 
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TABLE  A- 2  DIMENSIONLESS  DATA  (Continued) 


BRASH  ICE 

R  v  R  v 2 


41 

1 

.462 

.488 

.491 

.238 

42 

1 

.449 

.734 

.478 

.  .39 

43 

.413 

.446 

.439 

.199 

44 

1 

.506 

.825 

.538 

.681 

45 

Qun 

i 

.416 

.126 

.442 

.016 

46 

.311 

.445 

.331 

.198 

47 

Q0  * 

.172 

.470 

.183 

.221 

48 

.381 

.720 

.405 

.518 

♦Deleted  from  analysis,  see  Table  A-l. 
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TABLE  A-3  DATA  IN  FULL-SCALE  UNITS  (Continued) 
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♦Deleted  from  analysis,  see  Table  A-l. 


TABLE  A-3  DATA  IN  FULL-SCALE  UNITS 
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APPENDIX  B 
FRICTION  DATA 

TABLE  B.l  FRICTION  AND  ROUGHNESS  MEASUREMENTS 


TABLE  B.1  FRICTION  AND  ROUGHNESS  MEASUREMENTS 
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TABLE  C. 


APPENDIX  C 
ICE  DATA 

1  DAILY  ICE  PROPERTIES 


TABLE  C.1  DAILY  ICE  PROPERTIES  (CON'T) 


APPENDIX  D 


REGRESSION  ANALYSIS  OF  DIMENSIONLESS  DATA 


TABLE 

SOLID  ICE 

TABLE 

BRASH  ICE 

D.l 

«o 

D.8 

«o 

D.2 

D.9 

D.3 

oi 

D.  10 

QU 

D.4 

D.  11 

D.5 

^2 

D.  12 

§2 

D.6 

«3 

D.13 

«3 

D.7 

«4 

D.14 

«4 

D.  15 

1975  Data,  1.5  Ft. 
Solid  Ice,  Q0 , 
Model  2 

Definition  of  Variables 


0 

Solid  Ice  Intercept 

Brash  Intercept 

1 

R  .  »• 

R 

WBhz 

Pj?Bhz 

2 

-2-  -  y' 

V 

Tgh 

Sgh 

3 

°  =  o' 

y2 

gh 

4  y1  x  o' 

The  computer  printouts  contained  in  this  Appendix  are  the 
result  of  forced  multiple  linear  regression  to  the  three  pre-selected 
predictor  equation  forms. 


D-l 


TABLE  D.l 


RE  SI  S  T  Ft  N  C  E  FOR  0  0 


MULTI RLE  LI NEAR  REGRESS  I DM 


NUMBER  Or  OBSERVATIONS  16 
NUMBER  OF  VARIABLES  4 


INDEX 

1 

MEANS 

1 1 . 3958 

STD  DEVS 
3.84448 

C 

J 

4 

181. 6  5‘  6 
116.718 

96. 1816 

1  0 1'  ■  -■  1  ci 

CORRELATION  COEFFICIENTS 

Ql.i  NUMBER 

1 

1 

.674056 

. 638834 

. 846 1 43 

0 

ROW  NUMBER 

. 674056 

a 

i 

. 1 08569 

.  7 46696 

0 

ROW  NUMBER 

. 638834 

.10857 

1 

.644763 

0 

ROW  NUMBER 

.848143 

4 

. 748696 

. 644783 

1 

0 

SUNS  OF  SQUARE' 

S  AND  CROSS 

PRODUCTS 

ROW  NUMBER 
1780. 96 

1 

9b. 1833 

87438.5 

19774.3 

0 

ROW  NUMBER 

96. 1338 

d 

6 . 88369 

1400.44 

1809. 13 

o 

ROW  NUMBER 

67 43  5 .  ”> 

3 

1400. 44 

•  49763 IE  06 

. 38691 4E  06 

0 

ROW  NUMBER 
19774.3 

4 

1809. 18 

. 3869 1 4E  06 

.  86  78 06 E  u6 

0 

kt-M*  A'il»  I 

D-2 


1 


TABLE  D.2 


RESISTANCE  C0R  C'4B 
MULTIPLE  LINEAR  PEERESS I Uh 


NUMBER  QF  OBShRVHT  I  DM.: 
NUMBER  □  r  '••‘RRIftBLES 


INDEX 

1 


ME  RM  S 
1 1 . 1 EE 
.7176 
1 06. 78 
74.E516 


STB  DEVS 
1.43713 
. 33 1 333 
E 1.73 13 
41 . 64* 


!.!jR:R£LRTIDn  COEFFICIENTS 


Ruiii  NUMBER 

1 

pQi.i.l  NUMBER 

. 365653 

ROW  NUMBER 

-. 163357 


. 365*53 


1 


-. 16335 


-.355136 


-.355136  1 

. 0534  . 33EE44E- 0 1 

SUM S  OF  SG'UPR'ES  RND  CROSS  PRODUCT  .> 


ROi.,t  NUMBER 

•  *3*3 '5*  2  St* 


ROW  NUMBER 
6E7 .301 

ROW  NUMBER 
43.  04 

ROW  NUMBER 
5313. 36 

ROW  NUMBER 
4355 . 43 


1 


43 .  04 


3 1 8 . 36 


158 IE  371.303 


t-  -*  \  •*;  f  i  ~j,  58303 

•-•I  X  •  L  -  - 


334.381  33747.8 


0-4 


4  355. 


344 


34  5  04 


i 


i 


330534 

:44E-Ol 


TABLE  0.2  (Continued) 

REGRESS IQM  N'JNEER  1  DEPENDENT  VARIABLE  IS 


INDEft 

B 

STD  ERROR 

0 

8. 7 1938 

.567901 

4 

. 324 13 OF- 01  .68 

■  6  33E-  08 

STANDARD  ERROR 

OF  ESTIMATE^ 

. 569411 

COEFFICIENT  OR 

VARIATION*  .511 

734E-01 

R-SQUHRED= 

.  3  8  EG' 6  R= 

.  9398:5 7 

ANALYSIS  OF 

VAR  I  AriCE 

DR 

PIS 

REG  1 

i'  ■  C.  *1'  o  £' 

7. 83864 

RES  3 

•  'r*  ?cL*z>  o  7 

.  384889 

ACTUAL 

PREDICTED 

RES  I DUAL 

1  0 . 3 1 

1 0.4591 

849058 

1 3 . 3 

11. 4c'S 

. 3 7 1 993 

9.  E9 

9. 48099 

-. 13099 

1 1 .  07 

1 1 . E3 33 

168164 

1 E .  7  6 

13. 0838 

-  .  I'd  3  r  AT 

T-RATIG 
.  5 . 3535 
4.7413 


DURE  I N-MAT  SON  S  T  A  T .  =  3 .  35357 

REGRESSION  NUMBER.  3  DEPENDENT  VARIABLE  IS  1 


INDEX 

0 


3.51331 
3. 63331 


STD  ERROR  T-RATIO 

.  4435  17  13.  '?G  (' 9 

. 564477  3.43736 


STANDARD  ERROR  OF  ES7IMATE=  .431141 

COEFFICIENT  OF  VARIATION*  .3875 USE- 01 


•  — _  nUAR'ED— 


938499  R  — 


D.  F.  = 


ANALYSIS  OF  VARIANCE  TABLE 
SS  MS  F 

7.  7  0  368  7.7  0353.  41. 4433 

.55734 3  .  1 8 5 6 3 3 


PREDICTED 


RESIDUAL 


t 

1 0.8513 

— . 64 1 38 

13.1 849 

.  1  15077 

- 

9. 15738 

. 138134 

l 7  35 

.  3  35  0  06 

18.  7  m  i_.  9 

. 59 068  0£- 0 1 

1 .33-: 


D-5 


TABLE  D.3 


R E S I  ST Hf  iC  E  F 0 P.  0 3 ft 
MULTIPLE  LINEAR  REGRESSION 

HUMBER  OF  OBSERVATIONS  5 

NUMBER  QF  VftRIfiBLES  4 

STD  DEVS 
3. 34933 
. 430099 
41.6031 
5  7 .011 8 

SOP  PELft V I OH  CDEFP I C I  EH  TS 


ROW  NUMBER 

1 

■t 

1 

. 63  7 453  . 7 447 28 

. 9 084 06 

o 

R  u  10  H  b  M  B  E  R 

5 

. 337453 

1  .£~8631E-Ul 

.  oo  9  8  0  8 

0 

P1  Qi.it  NUMBER 

. 744733 

. 273633E-01  1 

. 482379 

0 

ROM  NUMBER 

4 

. 903403 

. 8 3 9 9 0 9  • 438879 

1 

0 

SUMS  QF  SQUARE 

S  AND  CROSS  PRODUCT'S 

RDM  NUMBER 

1 

507. 9 

32.9145  5230.33 

36  04 . 5>  2 

0 

F.Oi.0  NUMBER 

32. 9145 

2.55  07 6  3 02 .534 

234 . 42 

o 

RQi.O  NUMBER 

5220.33 

8 02 . 534  “«8  3  0  i .  4 

34792.2 

o 

RQlit  NUMBER 

4 

36  04 . 52 

284.42  347 92 . 2 

31301 . 5 

Ij 

INDEX  MEANS 

1  9 . 622 

3  .6018 

3  99.33 

4  6  0 . 49-5(3 


TABLE  D.3  (Continued) 


REGIES SION  NUMBER 


DEPENDENT  VARIABLE  IS 


I ND  E.-. 
0 
4 


E  STD  ERR OF 

6 .  3929  1 .  1  88  1 7 

532752E-  01  .  141  327E-  0 1 


T-PATIO 
5. 69698 
3 .  7  E  S  -J  — 


standard  error  df  e  : T  i vHTt.= 

COEFFICIENT  OF  VARIATION* 

R-  SQUARED*  .  8258  07  R  = 


1 . 6 1 7 1  7 

ICS  <> " 

. 903409 


D .  F .  = 


hmhl  V  _  I  S  OR  ‘•■'ARTANuE  TABLE 


REG 

FES 


3.  .  0  99 
7.345b' 


m; 

7  .  I  .  -  -4 ‘-4 

c.bl  528 


14,  lb  Sc 


Fib  TL'HL 
5. 63 
9  5 


14 


PREDICTED 
7. 0240b 


95 


7 . 66494 
1 4 . 6333 


RE . I DUAL 
-  1 .  394  r>- 
-  .  1 6  5  4 b 
-.  1 1223  3 
2.23507 
c  7  7  E  '4  .  E  -  0 1 


DO  RE  I N— filFiTSGN  C  i  AT  .  = 


1 . 43 1 96 


REGRESS  ID N  N u M B  E  R 


DEPENDENT  VAR  I  ABLE  I! 


INDEX 

0 


B 

6.3998 
5 . 35488 


STD  ERRGP 
2. 3 3245 

•2*  ■  C.  T’  ►£> 


T-RATI1 
2-  •  i  4  :  'z  1 
1  .  b  9  9 b 


STANDARD  ERROR  OF  ESTIMATE* 
COEFFICIENT  OF  VARIATION* 


"-SQUARED* 


4  7  8 €•  0  2  R* 


2 . 8  09  0 

. c  9 1 —48 

.  68  7 4b  1  0 .  h  .  * 


ANALYSIS  OF  VARIANCE  TABLfc 


REb 

:c 


DF 

i 


n.z.-i 


d  1 . 2 1 3 
2  3.6725 


MS 

Cl  1  .  Cl  1  -• 

7.  8908' 


AC  TUAL 
5 . 63 


9  •  86 

9.  95 

14.  id 

DURBIN-MAT SGN  ST AT.* 


PREDICTED 
7 . 8  i  00.8 
9.7141 
11. 74. -.4 
r'  .  88*61 
1 2 . 0-39 


RESIDUAL 
- 1 . 68  0  08 
-1.7641 
-  1 . 98-33  7 
c .  6b 1 39 
2 . 666  08 


.  8  <  -''484 
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TABLE  D.4 


RES  I 

STANCE  FOR  91 

MULT I RLE 

LINEAR  REbFbSSIGN 

NUMBER  OF 

OBSERVATIONS 

Q 

number  of 

VARIABLES 

4 

I  NIi  EX 

MEANS 

STD  DEVS 

1 

i 

1  i.i  .  76 1  8 

8.977  03 

CL 

.  K>  3  06  39 

.400487 

3 

169. 389 

61.1 733 

4 

1 08 . 7  09 

8  0 .318 

C  OF  P 

ELATION  COEFFICIENTS 

ROM 

r^ur*i  E;c£P 

1 

i 

.913474 

. 883889  . 89984 

ROo.' 

NUMBER 

£ 

. 913474 

l 

83889 IE- 01  .35611 

s  □'ll 

NUMBER 

3 

. 8 8 38 89 

.  8 3 8 '90 EE -01 

1  .49491 

P  □  K.l 

NUMBER 

4 

.  89  r*8-4“ 

.856115 

.494916  1 

sum: 

OF  SOURCE: 

S  AND  CROSS  f'RODUUTS 

PQH 

NUMBER 

1 

1 1 14.41 

69.3511  16 

889  18856.1 

ROW 

NUMBER 

£* 

69 . 85 1 1 

4 . 86498 

978.83  837.503 

saw 

NUMBER 

165 

c'9 

978.831 

cod 1 76E  06  . 135 180 E  0 

r'  Ll  l»J 

NUMBER 

4 

1  cc 

56.  1 

337.508 

135180E  06  . 157 958 E  0 

TABLE  0.5 


RESISTANCE  FOP  Q3 


MUL  T I PLE  L I NEHP  R EGRESS I DM 


NUMBER  OF  OBSERVATIONS 


NUMBER  OF 

VAR  1 RBLES 

4 

INDEX 

1 

MEANS 

9.74375 

.53775 

STD  DEVS 

5. 17008 
.418738 

4 

137.313 

86. 0573 

73. 534c 

i’ll'. 
— i  ; 

RELATION  COEFFICIENTS 

ROW 

NUMBER 

1 

1 

. 937439 

. 550737 

. 954 1 96 

ROW 

NUMBER 
.  937439 

1 

.406931 

.965133 

ROW 

NUMBER 
•  5  *'  0  r  C  r*' 

. 406931 

1 

. 586  7 66 

ROW 

NUMBER 
.  954197 

4 

.965133 

. 536 7 63 

1 

SUM 

3  OF  SQUARES 

AND  CROSS  c' 

RDD'JC  TS 

ROM  NUMBER 
946 . 633 

ROW  NUMBER 
60.  0333 

ROW  NUMBER 
1 1430. 9 

ROW  NUMBER 
4-347 . 5 1 


11430.  9 


9347.51 


4 

613.665 


683 . 573  613. 665 

.159886E  06  . 1 05455E  06 

. 105455E  06  97097.7 


REGRESSION 

NUMBER  1  DEPENDENT  VARIABLE  IS 

1 

INDEX 

B 

STD  ERF-DR 

T-RATIO 

o 

3 . 97  03  3 

. 946073 

4. 19664 

4 

. 67083 IE —  01 

858747E— hi 

7.31233 

STANDARD  ERROR  OF  ESTIMATE* 

1  •  6  r  U  r’  c! 

COEFFICIENT 

OF  VARIATION* 

. 171466 

A— SQUARED* 

.310491  A* 

.954197 

D .  F .  = 

6 

ANALYSIS 

OF  VARIANCE 

TABLE 

DF 

ss 

MS 

F 

REG  1 

170.  36 

170.  36 

61 . 0325 

RES  6 

l^. 7473 

2.7913 

ACTUAL 

PREDICTED 

RESIDUAL 

5.6  3 

5. 36345 

-.838454 

17.2 

13.0979 

-.897949 

4.46 

3.5? 

5. 595 

9. 7919  7 

-1.135 

-1.20197 

4.  £S 

4 . 82658 

— . 54657 6 

8,  i.i? 

7. 15152 

•  9  i  3  4  o  «» 

1 4 .  3 
15.4 

10.8587 

15. 7599 

■  441  J  - 
-.359861 

DURBIN— ilATSDN  ST  AT .  = 


1 . 4231 


R E6RE 3 3 1 □ N  N U M £  £  R 


DEPENDENT  VARIABLE  IS  1 


INDEX 

B 

STD  ERROR 

T-RATIO 

0 

2.94055  . 

‘  1.839 

c  •  o  r  -■  c! 

C. 

1 1 . 575 

1.75419 

6 • 5984“ 

STANDARD  ERROR 

OF  ESTIMATE* 

1 . 94243 

COEFFICIENT  OF 

VARIATION* 

. 199454 

R-SQUAREB* 

. 878885  R= 

. 937439 

D .  F .  =  6 

ANAL VS I S 

OF  VARIANCE 

ThBLE 

DF 

SS 

MS 

F 

REG  1 

164.446 

164.446 

43. 5398 

RES  6 

22 . 66 1 5 

3. 77692 

ACTUAL 

PREDICTED 

RESIDUAL 

5. 63 

4.93145 

. 698549 

17.2 

16. 3212 

. 878765 

4.46 

■v .  89673 

— .  436 726 

8.59 

10.2558 

- 1 . 66594 

4 . 28 

4.7231 

-.443101 

3.  09 

9 . 596 1 7 

-1.50617 

/ 

14.3 

10.5222 

3.77733 

15.4 

16. 7032 

-1.30321 

CURB IN— WAT SON  STAT 


2.63165 


D-ll/D-12 


TABLE  D.6 


RESISTANCE  C0R  0  3 
MULTIPLE  LINEAR  REGRESSION 


NUMBER  OR  OBSERVATIONS 
NUMBER  OF  VARIABLES 

INDEX  MEANS 

1  9.545 

£  . 655£5 

3  159.575 

4  100.769 

CORRELAT I ON  C OEFF I C I ENTS 


4 

4 

STD  DEVS 
2.80551 
. 372673 
16. £551 
5  1 . 75.-i>: 


RQI.it  NUMBER 

1 

Rou  number 

. 596 1 o 

rou  number 

-.523373 

ROM  NUMBER 


.  '398 1 6 


— . 3 £6  073 


•  ••  —  u_  . 


d— c*  U  r  r' 


.983056  -.  7: 

SUM 3  0-  SQUARES  AND  CROSS  PRODUCTS 


9£6£35 


— .  ('  568 


R  0  U  N 1  j  M  B  E  F 
388. 054 

ROM  DUMBER 
5!  {  .  6  554 

Rum  NUMBER 
6080.  35 

ROM  NUMBER 
4850. 93 


1 


ci.  8-554  6  0£  0 . 95 


£.13407  403. £33 


4£5 0 . 93 


3£1 . £86 


4  05 


-•  •  -J 


•  1  Oct  49E  0*6  6 £4-i 6 .  -5 


3£ 1 . £36  68456 . 3 


4365£.  9 


D-13 


i  \ 


TABLE  D.6  (Continued) 


REGRESSION  NUMBER  1  DEPENDENT  VARIABLE  IS  1 


I N  D  £  7. 

B 

STD  ERR DR 

T-RATIO 

0 

4. 48394 

1 . 59393 

£.31 3  04 

4 

. 30££44E— 01 

.  14453  IE- 01 

3.475 

i TAN  DA r d  e r r 0 P 

□F  EST IUATE= 

1 . £9557 

COEFFICIENT  OF 

VARIATIONS 

. 135733 

R  -  S  £.'  U  AR£D= 

.85791  R= 

. 9£6£ 35 

D .  R .  =  £ 

ANALYS  I S 

QC  VARIANCE 

TABLE 

DR 

NS 

p 

REG  1 

£0. £691 

£0. £691 

12. 0756 

RES  £ 

3. 35703 

1.67851 

HCT URL 

predicted 

RESIDUAL 

6  .  39 

5 . 99  0  07 

. 399532 

1  0 . 64 

1 0. 0736 

. 566439 

H  .  -1 

9 . 88 0  75 

-1.58075 

15. 35 

1£. £356 

. 614382 

ItfjRB I N— WATSON  ST A 7.=  5.81? 


REGRESSION  NUMBER  £  DEPENDENT  VARIABLE  IS  1 


INDEX 

0 


B 

5.  11: 

6 .  76: 


STD  ERROR 
1 .70998 
£. 34109 


T-RATIG 
990E8 
!.  83897 


STANDARD  ERROR 

OF  ESTIMATES 

1.51115 

COEFFICIENT  CF 

VARIATIONS 

. 158318 

P-30UARED= 

T  1  1  f-i  f*i  y  1  = 

. 8 98 16  D . R . = 

ANALYS I S 

OR  VARIANCE  TABLE 

DR 

3  - 

MS 

REG  1 

19. 059 

19.  059  8.. 

RES  £ 

4.56714 

£ .  £'8357 

ACTUAL 

PREDICTED 

RES  I DUAL 

6. 39 

&.  £766£ 

. 1 1 333 

1  0.64 

9. 47568 

1 . 1 643£ 

■  2’ 

1 0 . 03 03 

- 1 . 7  3  027 

15. 33 

l£. 3974 

. 452567 

B'JRBIN-MHTSQN  SThT 


3.11 966 


TABLE  D.7 


n 


RESISTANCE  p0R  Q4 

multiple  linear  regression 


NUMBER  OF  OBSERVATION'! 
NUMBER  OF  VRRIRELES 

INDEX  MEANS 

1  10.478 

8  .6504 

3  £30. 56 

4  151.676 

CQRRELR T I ON  COEFF I C I EN T S 


STD  DEV; 
4.34373 
. 467844 
84. £484 
l S5. 847 


ROW  NUMBER 

1 

ROW  NUMBER 

.  33 1 073 


.881073 


1 


•GW  NUMBER  3 

. 4 8 3 8  0 6  . 5 133 0  0E-01 


.  48  38  06 


. 51 3300E-01 


1 


RGw  NUMBER 

. 364646 


4 


. 355 05  . 464  378 

SUMS  OF  SQUARES  AND  CROSS  PRODUCTS 


. 364646 


i  0  4 


.464378 


ROW  NUMBER 
646. 348 

ROW  NUMBER 
43.5346 

ROW  NUMBER 
18335 

ROW  NUMBER 
1 0343. 9 


48.5846  18385 


3. 06 v 07  c;Si 


1  034  3.  '3 


7  03 . 8  8 '3 


703.88 


. 884 1 8 1 E  06  . 1 9456 7E  06 

.  1 '3456  7b  Ob  .  1 7 8 8  r  8 1  06 


D— 1 5 


TABLE  D.7  (Continued) 


REGRESS lDh  HUMBER 


DEPENBEN T  VARIABLE  IS 


INDEX 

0 

4 


B  STD  ERROR  T-RATIO 

4 . 7232 7  1.13  0  3  7  4.1 7353 

,  3734  09E  —  U 1  .  5  9 'S 4jcE-  0 c  6- .  3  .■  3 So 


STANDARD  ERROR  OF  ESTIMATE® 
COEFFICIENT  OF  VAR I AT I ON® 


R  —  S  Q  U  HF'ED— 


9 S'  OS 4 4  R= 


1 . 5 Ob S3 
. 143756 

•964647  D .  F .  = 


REG 

W  “  ,\ 


DF 

1 


ACTUAL 
6 . 77 
1  3 .  0  7 
1 6 .  3 
4.  09 
i  c. .  lb 


ANALYSIS  '3C  VARIANCE  TABLE 
SS  MS  c 

91.1336  40.1933 


91.1936 
6. 30664 

PREDICTED 
6 . 5697 
11. 669 
17.5333 
5.54193 
11.0761 


RESIDUAL 
.  3  0  0  3 
1 . 40097 
- 1 . 33336 
-1.45193 
1. 03391 


DURB IN-MATSON  3 T Fi T .  = 

REGRESSION  NUMBER 


13301 


DEPENDENT  VARIABLE  IS 


1 


INDEX 

0 


B 

4 . 6  b  7  7 
3. 93947 


STD  ERROR 
3.17  0  0 1 

3 .  »■'  r’  Ob  !*" 
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APPENDIX  E 

METHOD  OF  SCALING  BUBBLER  SYSTEM  CHARACTERISTICS 


The  laws  of  similitude  governing  the  modeling  of  the  air 
bubbler  system  has  been  investigated  using  two  different  approaches. 
Each  approach  begins  with  the  fundamental  requirement  of  ship  model 
testing  that  forces  must  be  scaled  by  Xs  where  X  is  the  geometric 
scale  ratio;  that  is, 

F 

prototype  _ 

*model 

X  =  ^prototype 
£model 

where 

F  -  force 
L  -  length 

The  first  approach  was  developed  by  Lewis  [7].  He  assumes 
that  the  geometry  of  the  air  plume  is  similar  and  that  the  drag 
coefficient  CD  will  be  equal  in  model-scale  and  full-scale.  This 
leads  to  the  following  relationship: 


(E.l) 

(E.2) 


^prototype  _  ,5/2 
Snodel 

where  Q  3  volumetric  flow  rate  of  air. 


The  second  approach  is  based  on  the  work  of  Kobus  [8],  For 
air  discharging  through  an  orifice,  he  develops  the  following 
expression  for  the  buoyant  force  acting  between  the  orifice  and  a 
cross-section  at  x : 

0  -  f*)  (E.4) 


(E.5) 
E-l 


F(x)  »  -  In 

% 

h*  *  h  + 


where 


u b  =  mean  rising  speed  of  air  bubble  stream 
x  =  vertical  distance  above  orifice 


F  =  buoyant  force 

Q  =  volumetric  flow  rate  of  air  (STP) 

0 

Patm  =  atm0SP^eric  pressure 
»  mass  density  of  water 
g  =  acceleration  due  to  gravity 

Kobus  [8]  also  states  that  the  mean  rising  speed  of  the  bubble 
is  independent  of  orifice  size,  spacing,  and  geometry  and  is  a 
function  of  air  flow  rate  only  as  follows: 

ub=  K(Qo)0-15  (E.6) 

where  K  -  constant. 

The  ship  modeling  laws  require  that  this  buoyant  force  scale 
according  to  the  following  relationship: 


[ 

f(x)] 

prototype 

[**(*) 

]  model 

Substituting  equation  E.6  into  equation  E.4  and  forming  the 
ratio  in  equation  E.7  results  in  the  following: 

[(a,)*,,!  fa  0  -  xA«)]  prototvpe  _  ^ ,  q 

[(«.)••“  ('  '  *"■*>]  model 

Solving  for  the  ratio  of  air  flow  rates  leads  to: 


(^prototype  =  ?3.53  [*”  ^  x/h*^}  prototype  (E  q\ 

(Vmodel  [in  (1  -  xA*)]  ^ 


E-2 


Equation  E.9  contains  x  and  h  ,  and  these  cannot  be 
eliminated  because  P  .  cannot  be  scaled  for  the  model  tests. 

atm 

A  comparison  of  the  scaling  of  the  air  flow  rates  by  equation 
E.9  and  equation  E.3  is  made  in  Table  E.l  for  the  l/24th  scale  model 
of  the  WYTM.  It  can  be  seen  that  in  this  case,  the  scaling  using 
equation  E.3  results  in  higher  full-scale  air  flow  rates  for  a 
measured  model -scale  air  flow  rate. 

Since  these  formulations  have  not  been  verified  with  full- 
scale  air  bubbler  systems  on  ships  in  ice,  good  engineering  judge¬ 
ment  dictates  that  the  scaling  be  done  according  to  the  relationship 
which  predicts  the  higher  air  flow  rates  for  the  full-scale  ship. 

It  is,  therefore,  recommended  that  equation  E-.3  be  used  as  the  scaling 
law  for  air  flow;  that  is: 


^prototype  =  xs/2 
Snodel 


(E.10) 
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TABLE  E.l 

RATIO  OF  AIR  FLOW  RATES  FOR  THE  140  FOOT  WYTM 


{q  ) 

v*o 'prototype 
model 


X 

Eq.  E.8 

Eq.  E.2 

(feet) 

2 

2399* 

2822 

6 

2280 

2822 

12 

2086 

2822 

'•**  h*patJW 

=  12.0  +  33.9 

=  45.9  ft. 
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